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We consider multidimensional gravity with a Lagrangian containing the Ricci tensor squared and the 
Kretschmann invariant. In a Kaluza-Klein approach with a single compact extra space of arbitrary di- 
mension, with the aid of a slow-change approximation (as compared with the Planck scale), we build a class 
of spatially flat cosmological models in which both the observed scale factor a(r) and the extra-dimensional 
one, 6(t) , grow exponentially at large times, but 6(r) grows slowly enough to admit variations of the ef- 
fective gravitational constant G within observational limits. Such models predict a drastic change in the 
physical laws of our Universe in the remote future due to further growth of the extra dimensions. 
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1 Introduction 



Multidimensional models of gravity and cosmology 
as low-energy limits of unified theories of physical 
interactions [1] are a powerful tool for studies of the 
present challenges to modern physics [2-4], such as 
the dark energy and dark matter problems, possible 
variations of fundamental constants [5-7], the role 
of strong field objects (black holes, wormholes) etc. 

It has been recently argued [8-10] that multidi- 
mensional gravity with curvature-nonlinear terms 
in the action can be a source of a great diversity of 
effective theories able to address a number of im- 
portant problems of modern astrophysics and cos- 
mology using a minimal set of postulates. Among 
such problems one can mention the essence of dark 
energy, early formation of supermassive black holes 
(which is a necessary stage in some scenarios of 
cosmic structure formation), and sufficient parti- 
cle production at the end of inflation. In this ap- 
proach, the particular value of the total space-time 
dimension D > 4 and the topological properties 
of space-time are supposed to be determined by 
quantum fluctuations and may vary from one space- 
time region to another, leading to drastically differ- 
ent universes. Different effective theories can take 
place even with fixed parameters of the original 
Lagrangian. It can be shown that this approach. 
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without need for fields other than gravity, is able to 
produce such different structures as inflationary (or 
simply accelerating) universes, brane worlds [10], 
black holes etc. The role of scalar fields is played 
by the metric components of extra dimensions. 

One of the challenging problems of modern 
physics and cosmology is that of possible time- 
, location-, and scale-dependent variations of the 
fundamental physical constants, in particular, of 
Newton's gravitational constant G. Variable effec- 
tive constants are a common feature of multidimen- 
sional cosmologies, where these constants depend 
on the properties of extra dimensions which can 
vary from one space-time point to another. This 
problem, among others, was addressed in [8]. A 
number of examples of cosmological models were 
built, where the size of the extra dimensions was 
stabilized at a minimum of the effective poten- 
tial. In such models, the constants could vary at 
earlier stages while the effective scalar field only 
approached this minimum, but take on stable val- 
ues together with the scalar when this minimum is 
achieved. 

In the present paper, we discuss another type 
of multidimensional cosmologies (their possibility 
was also mentioned in [8]) in which the extra- 
dimensional scale factor can grow indefinitely at 
large times but is yet small enough at present. Such 
models are of interest since they predict rather an 
exotic, though sad, fate of the Universe we live in: 
sooner or later, its physical laws must drastically 
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change due to the growth of the internal space. 
We will show that such models can in principle be 
viable since they simultaneously predict an acceler- 
ated (de Sitter-like) expansion and, under certain 
conditions on the input parameters of the theory, a 
small enough variation of the gravitational constant 
in agreement with observations. 

2 The multidimensional theory 
and its reduction 

We consider a {D = 4: + di) -dimensional manifold 
with the metric 

ds^ = g^^dxf'dx" + e'^f^^'^hatdx^'dx^ (1) 

where the extra-dimensional metric components bab 
are independent of x^ , the observable four space- 
time coordinates. 

The -dimensional Riemann tensor has the 
nonzero components 

^ pa — ^ pai 

R-\a = e-2^i^''^rf + <5"^,/3^/3^ (2) 

where capital Latin indices cover all D coordinates, 
the bar marks quantities obtained from g^j^ and hab 
taken separately, = 5^/5 and 5°-^cd = ^i5^d~^d^c- 
The nonzero components of the Ricci tensor and the 
scalar curvature are 

^"p. = Rp. + di B'^, 
Rl = e-^PR'^ + 5l[ai3 + di{dp)\ 
R = R[g\ + e-'^PR[b] + 2(iin/3 

where {8(3^ = P^(3^' , □ = VV^ is the d'Alembert 
operator while R[g] and R[b] are the Ricci scalars 
corresponding to g^i, and bab , respectively. Let us 
also present, using similar notations, the expres- 
sions for two more curvature invariants, the Ricci 
tensor squared and the Kretschmann scalar /C = 
R^'^'^^'Rabcd : 

RabR^"" = R^uR^" + 2d{R^,B'''' + e-'^'^RabR"' 
+ 2e-^^R[h][ap + di{d(3f] 
+ di[n/3 + di(a/3)2]2 (4) 

/C = lC[g] + AdiB^^Bf"" + e-^f^Klh] 

+ Ae-^^R[h]{df3f + 2di(di - l)[{dl3ff. (5) 



Suppose now that bab describes a compact di- 
dimensional space of nonzero constant curvature, 
i.e., a sphere (/c = 1) or a compact di -dimensional 
hyperbolic space [11] {k = —1) with a fixed cur- 
vature radius tq normalized to the -dimensional 
analogue of the Planck mass, i.e., rg = I/ttij^ 
(we use the natural units, with the speed of light c 
and Planck's constant h equal to unity). We have 

R"'\d = kmlS^^d^ 
rI = kml{d^-l)5l 
R[b] = kmldi{di-l) = Rb. (6) 

The scale factor b{x) = in ^ is thus kept di- 
mensionless; Rb has the meaning of a characteristic 
curvature scale of the extra dimensions. 

Consider, in the above geometry, a sufficiently 
general curvature-nonlinear theory of gravity with 
the action 

S = \^d'^ j V^d^'xiLg + L^), 
Lg = F{R) + ciR^^Rab + C2/C, (7) 

where F{R) is an arbitrary smooth function, ci 
and C2 are constants, L.^ is a matter Lagrangian 
and = |det(5MAf)|- 

The extra coordinates are easily integrated out, 
reducing the action to four dimensions: 

S = ^V[di] ml j y/^gd'^x e'^^^ [Lg + L^], (8) 

where "^g = \ det(5'^jy)| and V[di\ is the volume of a 
compact di -dimensional space of unit curvature. 

Eq. dSD describes a curvature-nonlinear theory 
with non-minimal coupling between the effective 
scalar field (3 and the curvature. Let us simplify 
it in the following way (putting, for convenience, 
m-B = 1): 

(a) Express everything in terms of 4D variables and 
I3(x) ; we have, in particular, 

R = Ri + <p + h, 

R4 = R[g], /i = 2din/3 + di(di + 1) (9/3)2, (g) 

where we have introduced the effective scalar field 

(l){x) = Rb e~2'3(^) = kdi {di - 1) e-2^(^) . (10) 

The sign of (p coincides with k = ±1, the sign of 
curvature in the di extra dimensions. 

(b) Suppose that all quantities are slowly varying, 
i.e., consider each derivative (including those in 
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the definition of R) as an expression containing a 
small parameter e; neglect all quantities of orders 
higher than 0{e^) (see [8,12]). 

(c) Perform a conformal mapping leading to the 
Einstein conformal frame, where the 4-curvature 
appears to be minimally coupled to the scalar 0. 

In the decomposition ([9]), both terms /i and 
i?4 are regarded small in our approach, which actu- 
ally means that all quantities, including the 4D 
curvature, are small as compared with the D- 
dimensional Planck scale. The only term which 
is not small is (f), and we can use a Taylor decom- 
position of the function F(R) = F{(j) + i?4 + /i) : 

F{R) = F{<l) + i?4 + /i) 

~F(0) + F'(,^).(i?4 + /i) + ..., (11) 

with F'{(l)) = dF/d(l). Substituting this, and the 
corresponding decompositions of the expressions @ 
and dS]), into Eq. ([8]), we obtain, up to O(e^), the 
following effective gravitational Lagrangian Lg in 
Eq. (HD: 

L 



g - F'{<l))Ri + F{ct>) + F'{<l))h + c,<t)- 
+ 2ci(pn(3 + 2(cidi + 2c2){dpf 



(12) 

with = ci/di + 2c2/[di{di — 1)]. 

The action ([8|) with (fT2]l is typical of a scalar- 
tensor theory (STT) of gravity in a Jordan frame. 
To study the dynamics of the system, it is helpful 
to pass on to the Einstein frame. Applying the 
conformal mapping 

g^u ^gt.u = \f{<P)\9^u, f{4>)= e'^'^F' (13) 

after a lengthy calculation, we obtain the action in 
the Einstein frame as 
1, 



S 



L 



-y[di 



1 



R^ + -KEmd(^) 



(/(sign F')L, 

^E(0)+i. 



Lm = (signF') 



1 

2^ 



F'{(p) 
JF" 



2 -^^mj 



2d 



+ -di{di+2) 



+ 



El 

4(ci + C2) 



-(sign F') 



;[F{cj))+C, 



(14) 

(15) 

(16) 
(17) 



'F'(</.)2 

where the tilde marks quantities obtained from or 
with g^i, ; the indices are raised and lowered with 
Qfiu', everywhere F = F{(j)) and F' = dF/dcj) 
expressed in terms of (/> using pU]) . 
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3 The cosmological model 

Depending on the choice of F(R) , the parameter ci 
and C2 and the matter Lagrangian in the action ([Tj) , 
the theory under consideration can lead to a great 
variety of cosmological models. Many of them were 
discussed in [8], mostly those related to minima of 
the effective potential (|17p at nonzero values of 4>. 
Such minima correspond to stationary states of the 
scalar (j) and consequently of the scale factor b = 
of the extra dimensions (see also [13]). If the min- 
imum value of the potential is positive, it can play 
the role of a cosmological constant that launches an 
accelerated expansion of the Universe. 

Here, we would like to pay attention to one more 
minimum of the potential VEin ) existing for generic 
choices of the function F{R) with F' > and lo- 
cated at the point = 0. The asymptotic 
corresponds to growing rather than stabilized extra 
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such an asymptotic growth at late times may still be 
of interest if the growth is sufficiently slow and the 
size b does not reach detectable values by now. Let 
us recall that the admissible range of such growth 
comprises as many as 16 orders of magnitudes if the 
D -dimensional Planck length l/rrio coincides with 
the 4D one, i.e. about 10~^^ cm. This estimate 
certainly changes if there is no such coincidence. 

Let us check whether it is possible to describe 
the modern state of the Universe by an asymptotic 
form of the solution for (/> — > as a spatially flat 
cosmology with the 4D Einstein-frame metric 



dsl 



dr 



\t)dx'^ 



(18) 



where is the Einstein-frame scale factor. At 
small </>, assuming a smooth function F[(j)) , we can 
use its Taylor decomposition 

F(0) = -2A + + C02 + ..., (19) 

where the form of the first two terms is chosen to 
obtain in ([7]) multidimensional Einstein gravity in 
the first approximation va. R; A is the initial cos- 
mological constant. For the kinetic and potential 
terms in the Lagrangian (jl4p we then have in the 
main approximation with respect to (p: 

Ko = ^[di{di + 2) + 8{ci + C2)]; 

Ve ~ Ae-'^i^ = A[di(di - l)]'^i/2|(/,|^i/2. (20) 

Assuming Lm = (i.e., restricting ourselves to 
vacuum cosmologies), we can write two indepen- 
dent components of the Einstein-scalar equations 
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for (3{t) and a^it) in the form 



Of/ 



Ad 



1 



4Ko 



(21) 



(22) 



It is hard to solve this set of equations exactly. 
Let us, however, notice that our equations are of 
the same type as those appearing in studies of in- 
flationary cosmologies, and the field /3 plays the 
role of an inflaton. Hence it seems possible to ap- 
ply the slow-rolling approximation frequently used 
there: we suppose 



« 3^/3, 



<. 3A e" 
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(23) 



and drop the corresponding terms in Eqs. (|22p and 
([2T]). Then we can express a^/aE from (p2]) and 
insert it to ([2T]) . getting 



2AKo 

which is easily integrated to give 
/2A 



4 



48Ko 



{t-to), 



where to is an integration constant which can be 
eliminated without loss of generality by changing 
the zero point of the coordinate t. For the scale 
factor we obtain 



with 



■p/t 



Of, OC 



tP 



p = 48Ko/dl 



(26) 



(27) 



Substituting the solution to the slow-rolling con- 
ditions ([23]) . we see that both of them hold if 3p » 
1 , or, in terms of the input parameters of the the- 
ory. 



3p = 36 



di + 2 8(ci + C2) 



dj 



> 1. 



(28) 



For ci + C2 > , we have 3p > 36 . This verifies a 
sufficiently good precision of our solution. 
Passing over to the Jordan frame with 



dsl 



dr'^ - a^{T)dx^ 



[dt^ - aE^{t)dx'^] (29) 



(so that r is the cosmic time in the Jordan frame), 
due to -F'(O) PS 1 , we can put simply / = e^^^ k, , 
to obtain 



(30) 



6(r) 



2A 



/2Ar/p 



p 



where we have fixed an integration constant by 
choosing the zero point of the time variable r . 

A further interpretation of the results depends 
on which conformal frame is regarded physical (ob- 
servational) [14, 15], and this in turn depends on 
the manner in which fermions appear in the (so far 
unknown) underlying unification theory involving 
all interactions. We here restrict ourselves to the 
simplest and maybe the most natural assumption, 
that the observational frame coincides with the fun- 
damental (Jordan) one, in which the initial theory 
([7j) has been formulated. 

Then an immediate observation is that the ex- 
ternal scale factor a(r) grows exponentially, in a de 
Sitter manner, which conforms to modern observa- 
tions if one properly chooses the constants, namely. 



(24) 



(25) {p - l)V2A/p = Hq^ 2.3x 10 



-18 g-l 



7.25x10"" yr"\ 



(31) 



where is the modern value of the Hubble pa- 
rameter. 

The internal scale factor 6(t) grows much slower 
for sufficiently large di , while the volume factor h'^^ 
behaves like q^^^^'^/p _ xhe effective gravitational 
constant is known to change inversely proportion- 
ally to the volume factor, so that 

(32) 



G/G = -2V2A/p, 



(here and henceforth the dot means d/dr). The 
dimensionless parameter of G variation is 



5 = G/{GH^) = -2/{p-l). 



(33) 



The tightest experimental constraint on G vari- 
ation has been obtained, to our knowledge, from 
lunar laser ranging (LLR) data [16], namely. 



G/G= (2±7)xl0 



-13 



yr 



(34) 



Thus a viable cosmology should predict 5 < 10"^ , 
which according to ([33]) constrains our model pa- 
rameters by 



p=l2 



' di + 2 
di 



+ 



(ci + C2) 
dl 



> 100, 



(35) 
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It follows that our initial field model ([7]) with 
only F{R) is unable to satisfy the constraint (j35p : 
to do so, it is necessary to invoke the Ricci tensor 
squared or /and the Kretschmann scalar, with the 
input constants ci and C2 such that 

ci + C2>dl (36) 

We conclude that, under the condition ([36j) . our 
model with the asymptotic behavior (|3U|) is poten- 
tially viable since it combines a de Sitter-like expan- 
sion of the observable Universe with a slow enough 
variation of the effective gravitational constant. 

Three more observations can be added. First, 
comparing the constraints ([35]) and ([28|) , we see that 
our approximation works so much the better, the 
smaller is variation of G. 

Second, according to ([3T]) . the constant A is de- 
termined by the modern value of the Hubble param- 
eter and is thus approximately the same as the ob- 
served cosmological constant. So this model suffers 
the same fine-tuning problem for the cosmological 
constant value as the Standard model. 

Third, the size 6(r) of the extra dimensions, 
given by (j30p . remains uncertain since it depends 
on the arbitrarily chosen values of the time vari- 
able r (note that observable quantities depend on 
Hq rather than r). Thus one can easily satisfy the 
conditions that 6(t) must be much greater than 
the Planck length but be still below the observa- 
tional threshold (about 10~^^ cm). We also see that 
even an exponential growth of the extra dimensions 
can be slow enough to conform to the observational 
bounds on the stability of fundamental constants. 
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